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Abstract. We introduce non-linear diffusion in a classical diffusion advection model with non local 
aggregative coupling on the circle, that exhibits a transition from an uncoherent state to a coherent one 
when the coupling strength is increased. We show first that all solutions of the equation converge to 
the set of equilibria, second that the set of equilibria undergoes a bifurcation representing the transition 
to coherence when the coupling strength is increased. These two properties are similar to the situation 
with linear diffusion. Nevertheless nonlinear diffusion alters the transition scenari, which are different 
when the diffusion is sub-quadratic and when the diffusion is super-quadratic. When the diffusion is 
super-quadratic, it results in a multistability region that preceeds the pitchfork bifurcation at which 
the uncoherent equilibrium looses stability. When the diffusion is quadratic the pitchfork bifurcation at 
the onset of coherence is infinitely degenerate and a disk of equilibria exist for the critical value of the 
coupling strength. Another impact of nonlinear diffusion is that coherent equilibria become localized 
when advection is strong enough, a phenomenon that is preculded when the diffusion is linear. 



1. Introduction 

We consider the equation 

dtu = dl{u^) + de {uJ * u) , t > 0, 6* e [0, 27r], 

m(O,0) =UQ{e) 0e[O,2^], 

ti(t,0) =ti(t,27r) t>0, 

deu{t, 0) = deu{t, 27r) t > 0, 



(1.1) 

= deu{t, 27r) t > 0, 

where m > and J*u{t,0) = K f^^ sm{9 — ip)u{t,(p) dip, with uq even, no(-) > and f^^uo{0)d9 = land 

/■T > is a constant. In the following we denote xi = xi{u) = Jq^ u{0) cos{0)dO for any u € -^^"'^([0, 27r]). 
The right hand side of this equation is comprised of two components, with the first representing a 
nonlinear diffusion and the second one a nonlocal advection term. The dynamics of the solutions result 
from competition between these two terms, as the first one tends to spread solutions whereas the second 
one, on the contrary tends to concentrate them. In this work we examine the changes in the dynamics 
of the above equation depending on the parameters m and K, which control, respectively the diffusion 
and the strength of advection. 

For the special case of linear diffusion, i.e. m = 1, this equation arises in various contexts, including 
as a mean- field spin X-Y models |35|, a Doi-Onsager or Smoluchowski model for nematic polymers 
[15| [T6l I43| . or a Kuramoto or Sakaguchi model of synchronization [H [36l [2^. Its equilibria set and 
dynamics have been analyzed in great detail |28l[l6l[3Sl|7l|20l|371[I71[I][l5l[71ll3l|20^ The picture that 
emerges from these studies is that there are two distinct regimes depending on the value of K. On the 
one hand, for K < ^, all solutions tend to the trivial equilibrium u = l/(27r). On the other hand, for 
K > ^, the equation admits a unique non trivial equilibrium (up to a rotation) that attracts all solutions 
except those that lie on the stable manifold of the trivial equilibrium. The non trivial equilibria are 
called coherent because they have a single maximum on [— 7r,7r], corresponding to a region of maximal 
density, where the population described by u aggregates. In this sense, the bifurcation taking place at 
K = ^ separates a regime where diffusion dominates from the one where advection promotes coherence. 
This transition whereby the number of equilibria of the system changes as K crosses a critical value is 
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the key property of the model. This sudden change of dynamics accounts for phenomena such as onset 
of synchrony in coupled oscillators, and the transition from isotropic to nematic phase in Doi-Onsager 
models. We refer to it as a coherence transition. Its characterization is one of the main motivations for 
the large number of studies devoted to this model. In the present work, our purpose is to analyze the 
impact of nonlinear diffusion (i.e. m > 1) on this transition. 

Nonlinear diffusion equations, with or without advection, have been studied as mathematical models 
for many important phenomena, such as diffusions in porous media or spread of biological populations, 
that show original and complex dynamics [6l[2l[l0l[ig|32l[22l[23l[2il[5l[l8l[Il]. In the context of 
biological aggregation, equations such as (|l.ip . with space variable 6 G R'^, m > 1 and various coupling 
functions, have been subject to an intense study, including mathematical results on derivation and well- 
posedness of the equation ( |3H [271 E El HSl [30l [l2]), and several phenomena that were known for the 
Porous Medium Equation (PME) [2l|3l HI 02] such as finite propagation speed of interfaces and existence 
of travelling waves or stationary solutions [111 W\\ [3T1 [29l I13| . Degenerate diffusion (m > 1) has been 
introduced as a biologically realistic mechanism for the apparition of clumps in aggregation models, i.e. 
the apparition of small groups of individuals with sharp edges, which had only been previously observed 
in one dimensional models with highly specific coupling kernels. In case m = 2 and for a general class 
of coupling functions J, Burger and al. [Tl] have studied the existence or non existence of equilibria 
in any dimension [9 S R*^), and uniqueness, monotony and compact support properties of equilibria 
in dimension one {9 G R). Taking m = 3 and J{9) = Ke~^^^, and considering equation (ll.ip with 
£ R, Topaz and al. [41j have shown in particular that depending on the coupling strength, various 
equilibria or periodic solutions with compact support exist. Degenerate diffusion has been introduced 
and considered in several other diffusion advection reaction equations. See |391 \38\ |40] for the classical 
Keller-Segel model of chemotaxis with degenerate diffusion for example. Despite this large number of 
studies, to our knowledge, the effect of nonlinear diffusion m > 1 in equations of the form (II. ip has 
never been investigated so far. This is the aim for our paper. 

Our main results are summarized in the following theorem: 

Theorem 1.1. For all m > 1 and K > 0, equation is the gradient flow of a free energy and all 

smooth solutions converge to the set of equilibria of lil.l]) . Denoting K = except for 

m = 2 and K = 2 this set of equilibria is finite. 
For any m > 1 we have: 

• For all K < 2, the uncoherent equilibrium ^ is locally stable. 

• For all K > 2, there is a unique pair of (locally stable) coherent equilibria of whose bassins 
of attraction contain an open and dense set of initial conditions. 

When 1 < m < 2, for all K < 2 the uncoherent equilibrium is globally stable. 
When m > 2, there is a < Kc{m) < 2 such that 

• when K < Kc the uncoherent equilibrium is globally stable, 

• when Kc < K < 2 the dynamics of U.l\) is bistable: the uncoherent equilibrium ^ is locally 
stable, there is unique pair of (locally stable) coherent equilibria of and all solutions of 
U.l\) converge to one of those. 

When m= 1, the coherent equilibria of U.l\) are positive on [— 7r,7r] for all K > 0. However, for all 
m > 1, there is a Ki(m) > Kc{m) such that for all K > Ki, the coherent equilibria of are localized 
in [—IT, vr]. 

This theorem shows that nonlinear diffusion, besides localizing equilibria, modifies the scenario of 
transition to coherence. This modification is captured by the two-parameter bifurcation diagrams of 
equation (jl.ip in the plane {K, m) represented in fig. [1] (wherein the left panel is a magnification of 
a section of the right panel). These diagrams show that {K,m) = (2,2) where a highly degenerate 
pitchfork bifurcation takes place is an organizing center of the dynamics. Three qualitatively distinct 



3 



regimes labelled as (U), (B) and (C) come to meet at this point. These regimes represent, respectively, 
(i) the parameter range for which the uncoherent equilibrium ^ is globally asymptotically stable, (ii) 
the multistablity regime where ^ is locally stable and coexists with two pairs of coherent equilibria, 
with one pair being unstable and the other stable, and finally (iii) the regime of coherence where ^ is 
unstable and most solutions converge to either one of a pair of symmetrical stable coherent equilibria. 
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Figure 1. Equilibria bifurcation diagram. Solutions of converge to the unique 
(uncoherent) equilibrium ^ in region {U). In region (C), solutions typically converge to the 
unique coherent equilibrium, with support strictly included in [0,27r] in region iCi). In region 
(i?) one coherent and one uncoherent (locally) stable equilibrium coexist, m and K are the 

parameters of (fLTj) . and K = (^)~''™~^^- 



This paper is organized as follows. In section[2l we show that, as in the case m = 1 or in the PME case, 
equation (jl.ip preserves positivity and mass. In section[3]we introduce the free energy of (jl.ip . which 
allows to write equation (jl.ip as a gradient flow in the formalism of |34| |33] . A consequence is that all 
solutions of (jl.ip converge to its equilibria set when t — )• +oo. We further show that at each equilibrium 
point the linearization of (jl.ip is symmetric for an appropriately chosen scalar product. In section HI we 
give existence and uniqueness results for equilibria of Ol.ip and give analytical formula. This shows that 
the pitchfork bifurcation at ^ is supercritical when 1 < m < 2 (section l5.ip and subcritical when m > 2 
(section I5.2p . When m = 2, the pitchfork bifurcation is infinitely degenerate, we have K — 2 = o(x") 
for all n > 2 in a neighborhood of the bifurcation, where xi is the order parameter that characterizes 
equilibria in bifurcation diagrams, xi = corresponding to the uncoherent equilibrium and xi > 
to coherent equilibria. This bifurcation scenario and the transition from supercritical to subcritical 
pitchfork at m = 2 are shown in section [5^31 The bifurcation scenario in the limit m — )• +oo is discussed 
in section 15.41 Several remarks and discussions are in section [6j Details on numerical methods used to 
make figures of this text can be found in the appendix. 



2. Positivity and mass conservation 

The problem (jl.ip is reflexion invariant: if u{t, 9) is a solution of (jl.ip then v{t, 9) = u{t, —9) is also a 
solution of (jl.ip with vq{9) = uq{—9). In particular, if uq is an even function then any solution u{t, •) of 
(jl.ip is even. Here we focus only on even initial data and we suppose uq{9) = uq{—9) in the following. 
If Uq E L^i^) with Mo ^ 0, and u{t, ■) is an associated solution of (jl.ip on [0, T], then we have u{t, •) — 
and J^u{t,9)d9 = JgUo{9)d9 for all t G [0,T] (see lemma [1]). In the following we assume uq > and 
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Jgi Uo{9)d9 = 1. We use the notation 

^) = T7- + - ^nit) cos{ne), 

n>l 

and in particular we denote xi = xi (n) = u{e)cos{6)de. 

Lemma 1. Assume that uq € and u is a weak solution of jll-l]) in C^{[0,T], L^{[0,2tt])). Then we 
have jZ^u{t,e)de = JZ^uoi0)d0 for all t G [0,r]. Suppose that u £ C^QO, T], L^([0, 27r])) is a weak 
solution of with n(0, ■) = uq > 0. Then we have u{t, •) > for all times t G [0,T]. 

Proof For any time t e [0,T] and (/) £ C°°{[0,T] x [0,2tt]), we have 

/ u{t,e)(j){e)de = / uo{e)(i){e)de + / / {u{s,e)"'d'i(j){s,e) - u{s,e)j *u{s,e)d0(j){s,e) deds, 

Jo Jo Jo Jo 

and for (j){s, 9) = 1 we find /^^ u{t, 9)d9 = f^^ uo{9)d9 for aU t < T. 

Consider u^{t,9) = max(-n(t, 6*), 0) > 0. We have ii_ E Ci([0, T], ^^([0, 27r])) and u_ is a weak 
solution of (jl.ip . Then the first part of the lemma gives f^"" u.{t,9)d9 = J^"" u-{0,9)d9 = for all 
t e [0, T]. With n_ > this implies u(t, •) = and then u(t, •) > for all i e [0, T]. □ 



3. Gradient Flows 

3.1. Energy decay and convergence to steady states. In all this section, we assume m > 1. 
Equation (jl.ip can be written as a gradient flow for the Wasserstein metric: taking 

^(n) = [ u'^{9)d9 - — H u{9)u{ip) cos{9 - (p)d9d(p, (3.1) 

m-l Js 2 7^5x5 

we have formally 

5J^ m 
6u m — 
and 

de (^udg ^^"^^ ^ = de ^muu"^~'^ddU + Ku j sin(- - Lp)u{{p)d<^ , (3.3) 
so that (|l.ip is equivalent to 

The functional is a strict Lyapunov functional for (|l.ip . More precisely, if u{t,-) is a smooth 
non-constant solution of (II. ip . we have 



-u™-^ - ^ / cos(e - ip)u{ip)dip, (3.2) 



9tJ^(u(t)) = L»J^(M).atu = / -^^.dtud9 = - I u{t,9) 



5u{9) 

= - y n {mu"'~^dgu + iiTxi sin(-)) d9 < 
The solutions n € with deiu""'^) S L°° of 



d9 

(3.5) 



750(u"^-^)(0) + Kxi sm{9) = 0, (3.6) 

\?Tl — 1 / 

where xi = Jgu{9) cos{9)d9, are the equilibria of equation (jl.ip (see section 3] for a description of the 
set of solutions of (|3.6p ). 
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Proposition 3.1. The functional T : L™(S) R defined in [3J\) is . If u e L"^{§>) with u > and 
Jgud9 = 1, then we have 

K rn ,, K ^, . m ,, K , „, 



Proof. The regularity of is a direct consequence of its definition (j3.ip . The hypotheses it > and 
J = 1 imply in particular < xi{u) < 1, and the inequalities on J-'{u) follow. □ 

Theorem 3.2. Assume that u is a solution of with u{t, •) > and u{t, 9)d9 = 1 for all t > 0, 

and u £ C^([0,+oo[,L™(S)). Let C be the set of solutions of ( fT^) . 
Then we have 

distLoo(u(t),C) = inf \\u{t) — u\\ Loo — )• 0. (3-8) 

Remark. Lemma [T] and theorem 13.21 directly extend to initial data with uq > and f^^uo{6)d6 < 
+00 only. In particular, solutions of (|l.ip converge to its equilibria set even for non-even initial condi- 
tions. This also holds when m = 1, see [7J. 

Proof For ah t > we have \\u{t)\\'f^m < J='{u{t)) + § < J='{u{0)) + §, and we have G L°°{[0,+oo[, L"") 
or equivalently G L°°([0, -|-oo[, L^). By hypothesis and proposition 13. H the function t ^ J- o u{t) is 
C"'^([0, +oo[, M), not increasing and bounded below. Hence we have J='{u{t)) ^ I > -§ and 5* 
when t — t- +oo. 

Let t„ > be a diverging sequence and u„ = u{tn, •)• The sequence (un) is bounded in and (up 
to taking a subsequence), we can assume that n„ '^^^ (where denotes weak convergence in 
space). We also have 

/ Un (mu^~'^d0U + Kxi{un) sin(-))^ d6 — > 0, 
Js 

and mu™ dgu -\- Kxl{un)^/ynSin{■) is convergent in L^(S). We have Un — — > Uoo, so that xi(n„) — t- 
xi{uoo), and 

Kxi{un)\/u^sm{-) — > Kxi{uoo)\/u^sm.{-). 
In particular mu^ ^^"^doUn = ^ i dg (u^ ^^j is bounded in L^™ C L^. Then ^^'^ is bounded in 



2 

and compact in Since ^^"^ — >u^ ^^'^ weakly in L'"~i/2, we have ^^"^ — >u^ ^^"^ in 

^dgu^-^ +Kxism{9) 



and we deduce Un — >Uoo in The limit Uoo G satisfies f?e lioo [^ ^"^i ^^^oo ~'~ -^^i sin(0) j = 
at least in the sense of distributions. Since dgu^~^ £ L°°, we have Uoo € C. 

We have shown that the limit Uoo of any converging subsequence u{tn,-) is in the equilibria set C, 
and (ESI) follows. □ 



3.2. Linearization at coherent equilibria. Recall that we consider equation (jl.ip in a functions 
space with u > 0, J^ud9 = 1 and u is even: u{9) = u{—9). In this section, u denotes an equilibrium of 
equation (II. ip . 

Consider v, w such that Jg vd9 = Jg wd9 = 0. Suppose that they are V and W such that 

7; = a0(na0y) and w = dg{udgW). (3.9) 
Then we have (see \34\ |33] for example) 

/" = /" Ftyde = - f udeVdeWd9. 

Js Js Js 

For such V and we define the bilinear form (v, w) = — vWd9 = — /g Vwd9 (note that ( , ) depends 
on u). For such a u ^ we have {v,v) > 0, and then v,w iv,w) defines a scalar product. 
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Proposition 3.3. Let be the operator associated to the linearization of at the equilibrium u, 

LuV = dg (^mu^^^dgv + vJ * u + uJ * v) . 
Let V, w be even and chosen as above. Then we have 

{LuV, w) = -m f u"'-^{e)v{e)w{e)de +1.1* v{e)w{e)de = {v, l^w), (3.10) 



where J (9) = K cos{9). In particular, the operator L„ is symmetric for the scalar product 
Proof. We have 

{LuV,w) = / [mu"^~^dgv + vJ * u + uJ * v) deWd9. 
With de(u"^) + uJ *u = 0, this gives 

(LuV, w) = j {mu'^-^de {^i!^ +J*i^ udeWd9 

= -m [ u"'-^vwd9-m [ de{u"''^)vdeWd9 - [ J*vwd9. 



(3.11) 



Since ^, v, and OqW are even, we have J^do^u^ ^)vd0Wd9 = and the first equahty of (|3.1(jp 
foUows. We have 



J * v{9)w{9)d9 = JJ Kcos{9 - (p)v{(p)w{9)d^pd9 , 

hence the second term of (|3.10p is symmetric in v and w, and we have {LuV,w) = {v,Luw). □ 

The importance of the above result hes in that self-adjoint operators with compact resolvent have 
real pure point countable spectra, and they are diagonalizable in a Hilbert eigenbasis. The existence 
of spectral gap is then sufficent to show hyperbolicity and linear stability of equilibria, and classical 
stable, central and unstable manifolds existence theorem are available. This allows for the definition of 
an approriate projections and the study of small perturbations effects on the dynamics in neighborhoods 
of the equilibria set (see [3 1201 EI] i^i case m = 1). 



4. Structure of the set of equilibria 

To serve as a basis for comparison, we recall briefly the bifurcation scenario in case m = 1 (see 
|28| \TU[ [71 [2D]). When the coupling strength is below a critical value Kc, equation (|l.ip has a 
unique (stable) equilibrium ^ that is called uncoherent. When K = Kc a supercritical pitchfork 
bifurcation occurs, and for all K > Kc there are exactly three equilibria: and u{9) = l.^'^^'^^W ^ 
where c = J^^ e°"^°^^^^d9 and a = J^^ u{9) cos{9)d9, and u{- + vr). When K > Kc, this equilibrium u is 
locally stable and attracts all solutions of (II. ip such that J^^uq{9) cos{9)d9 > 0, whereas all solutions 
with J^^uo{9) cos{9)d9 < converge to u{- + vr). These two equilibria are non constant, they have a 
unique maximum on [— vr, vr] and they are called coherent in that sense. In particular, coherent equilibria 

have been shown to diverge from ^ as {K — Kc)^ (]T\\37\ [T7] and references therein). 

In this section, we establish the main elements of proof of theorem ll.il In the next section [5] we detail 
further bifurcations scenario (depending on m), which lead to the phase transition diagram of figure [T] 

4.1. The uncoherent equilibrium A pitchfork bifurcation occurs at the uncoherent equilibrium 
when K = Kf, = 2in (^)™ ^ ■ The uncoherent equilibrium is linearly stable when K < Kh and unstable 
when K > Ki^. In the following we rescale the parameter K = K^^ (^-^^ ^™ so that the bifurcation 
occurs a,t K = 2 independently of m. 
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Proposition 4.1. For all m > and K > 0, u{6) = ^ is an equilibrium of equation U.l]) . The 
linearization of equation at ^ is: 

/ I \ j{ 

dtv = Lj_v = m { — ] dgv + —xicos{9), (4.1) 

the operator Lj_ is orthogonal is the classical Fourier basis of i^(S) and its spectrum is given by 

27r 

Ai = ^ ~ 2m (^)™ and A„ = — n^2m (^)'^ for all n > 2. 
In particular, ^ is linearly stable if and only if 

K <Kh = 2ml — \ . (4.2) 

For K = Kh the operator Lj_ has a null eigenvalue, and for K > Ki, the equilibrium ^ is linearly 
unstable. 

4.2. Coherent equilibria of equation (jl.ip . 

Proposition 4.2. Suppose m > 1. Equilibria u of equation satisfy either u{0) = or 



m — 1 



u"-^ = Kxi[cos(9) + c],, (4.3) 

m ^ 

with J^^u{6)d9 = 1 and J^^u{9) cos{6)d9 = xi and c > —1, where x+ denotes max(a;,0). 

For every c > —1, there is a unique K > and xi €]0, 1] such that ( (^.g[ ) holds, they are given by 

ml 1 Jm{c) 

K = ^ , , ^ , , K and xi = ^ , , , (4.4) 

m - IJ^(c) /™(c)"^-2 I^(c)' ^ ' 

where Im{c) = f^^ [cos{9) + c]^ d9 and Jm{c) = jZ^ [cos{9) + c]^ cos{9)d9. 
In case < m < 1, the non-constant equilibria of equation are 

^m^l ^ l^^^.^ [c _ cos(0)] , , (4.5) 
m 

with c > 1, where xi and K are given by 

\Jm\y) , rn 1 ^ i \2— m 

X„^(c) l-mJm{c) 



with Zm{c) = Jl^ [c - cos(0)]|-™ d9 and J^{c) = /_^^ [c - cos{9)]]r"' cos{9)d9. 



Proof. Assume that m > 1 and n is a non-constant equilibrium of (jl.ip . Following equation (|3.6p . 
for all a £ [0, 27r] we have either u{a) = or u^'^^{9) = -^^^^-^xi (cos(0) + c) on an open interval 

containing a. With u'^~^ G W^'°°(S) and u > 0, this implies u'^~^{9) = —^^-^^—^xi (cos(0) + c)_,_ on 
[— arccos(— c), arccos(— c)], that is an increasing function of 9. Thus there is a unique c such that (j4.3p 
holds on Supp(ii) = [— arccos(— c), arccos(— c)] and u"^~^{9) = on [0, 27r]\Supp(u). 
We have 



m 



m — 1 



u{9)d9 = ( ^^"^ ^K i ] Im{c) = 1, (4.6) 



u{9) cos{9)d9 = ( — ^xi ) Jm{c) = XI, (4.7) 



m 

which give directly xi = j^. Then we have -fr^^ Jm(c)™'~^ = x™"^ and ([i^i]) foUows. 

In case < m < 1, equation (|4.5p is proved similarly. □ 
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Proposition 14.41 shows in particular that the branch of equihbria of proposition 14.21 appears by a 
pitchfork bifurcation at ^ when K = X^, or equivalently K = 2, and that it persists until K — )• +oo. 
Equation (j4.12p shows that when 1 < m < 2 the pitchfork bifurcation is supercritical, and when 
m > 2 the pitchfork bifurcation is subcritical. In both cases 1 < m < 2 and m > 2, the distance 

\l/2 

between coherent equilibria and the uncoherent equilibrium scales as I — 2 J in a neighborhood 
of the bifurcation. When m = 2 the pitchfork bifurcation degenerates, the distance between coherent 
equilibria and the uncoherent equilibrium is larger than some positive constant for any K > 2 (see also 
proposition 14.31 for the case m = 2 ). 

Proposition 4.3. In the case m = 2, (equations ^-4^ still hold) the equilibria are given by 

• for c>l, 

2 1 
K = — and xi = — , 

TT 2c 

• for — 1 < c < 1, 

2 arccosf— c) 

K = and xi = . 

arccos(-c) 2carccos(-c) + 2vl - 

Proof. In case m = 2, their respective definitions give /2(c) = 2carccos(— c) + 2\/ 1 — for |c| < 1 and 
/2(c) = ^ for c > 1, and J2(c) = arccos(— c) for |c| < 1 and J2(c) = vr for all c > 1. When m = 2 
equation (|4.7p reads 2K = Jm{c), and combined with xi = j^, this implies property 14.31 □ 
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Figure 2. as a function of c for m > 1. Coherent equilibria of (jl.ip are uniquely 
determined by K and xi, which depend on the parameter c g] — l;+oo[. When 1 < m < 2, 
Km{c) is decreasing on ] — l;+oo[, and for every K G [2;+oo[, there is exactly one coherent 
equilibrium of (jl.ip . This is illustrated by the top most curve on the panel (a) for m — 1.5. 
For m > 2, Km{c) has a minimum at some c* < 1, it is decreasing on ] — l,c*] and increasing 
on [c*]+oo[. The below most curve on panel (a) and the the four curves on panel (b) show 
that this holds for various values oi K — 3.0, K — 5.0, K = 8.0, and K = 15.0. For all m, 
we have K ^ 2 when c — > +00. When m > 2, for all K € [K„i{c*); 2] equation (jl.ip has two 
coherent equilibria. They appear by a fold bifurcation when K = Km{c*), and the unstable 
one disappears by subcritical pitchfork bifurcation at when K = 2 (see proposition 14.21 and 
illustrations in figure [21) 



Proposition 4.4. Assume that m > 1. 

• We have |xi| < 1 and K > for all c g] — 1; +oo[. 
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We have xi = 1 — + c) + o(l + c) and 



K 



l-m 



m + 1 



m + 1 



74-'"(l + c)-^+o(l + c)-^, 



(4i 



ly/ien c —7- —1, where the constants > and 7^ > are given in ^.16] ). In particular we 
have 



xi = 1 — fiK + o [K '"+1 I when K — )• +00, 



(4.9) 



W^th^l = Prn['^) ""+^7^ 

When c — )• +00, «;e /laue 

1 

Xl 



2 m— 1 
m+1 ^ 



1 



2(m - 1) c 
K = 2m (27r) 



(2 - m)(l - 2m) 1 

1 H TT, t; + O 



(m — 1)2 



l—m 



2-m 1 
(m — 1)2 



+ ^ 



and in particular 



K -2 = -{m - 2)xl + o {xl) . 



(4.10) 
(4.11) 

(4.12) 



Proof. We have Im{c) > for all c > —1, and since — )• (c + cos(0))^ is not increasing and non negative 
on [0, tt], we have Jm(c) > 0. Furthermore we have | Jm(c)| < Im{c) for all c > —1, and the inequalities 
l^il < 1 and K > Q follow. 
Denoting e = 1 + c, we have 



Im{c) 



arccos(— c) 



(c + cos(6'))™-i ^6* = 2 / (c + y)'"-i 



dy 



(4.13) 



2^e— / (1 - A) 



I m-l 



and similarly 



Jm{c) = 2y/ee 



1 

m— 1 



V^(2-eA) 



(1-A)' 



1-eA 



^A(2 - eA) 



dA. 



This implies xi — ^ 1 in (j4.8p , and with X 
precisely, we have 



ml 1 



m-l Xl 



we also have — t- +00 when c 



Xl = 1 — /3me + o(e) and K 



7i™e 2 (1 _ /3„e + o(e)) , 



m — 1 



where 



(1-A) 



1 

m — 1 . 



VA(2 - A) 



dX and /3 = - (1 - A) 



7 Jo 



I m-l 



VA(2 - A) 



zdA, 



(4.14) 
-1. More 

(4.15) 
(4.16) 



and this implies (|4.8p and (|4.9p . 



Using the Taylor expansion (1 + x)" = 1 + ax + ^^ x^ + °(" ^)(" 2) ^3 _|_ (2)(3;^), we find 



m— 1 



1 + - cos(6') de = 2ttc~ 



2-m 1 

1 + 77 Tv?^+0 

4(m - 1)2 c2 



when c — t- +00. Similarly, we have 



^11/ (3-2m)(2-m)l 
m — 1 c 



(4.17) 



(4.18) 



(m — 1)2 c2 

The estimate (|4.10p on xi is a consequence of (|4.4p and the two estimates above. The asymptotic of K 
is a consequence of this estimate on xi and K = (c))"'~i • Equation (|4.12p is a consequence of 

(ITOD and (liTT]l □ 
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The bifurcations in cases 1 < m < 2, m = 2, m > 2 and the transition at m = 2 are discussed in the 
following sections. The case m — )• +oo is also discussed in section [5741 

5. Global description of results 

After this analysis of the equilibria set of equation (jl.ip . we synthesize our findings by presenting, for 
various ranges of the nonlinear diffusion coefficient m, the corresponding bifurcation diagrams when the 
parameter K is varied. One of the main ingredients in our analysis of equilibria and their dependence 
on parameters has been the function c ^ K. It is the shape of this function that determines the number 
of non trivial equilibria of the system. The left panel in figure 2 illustrates how this function changes for 
various values of the parameter m. For 1 < m < 2, it is strictly decreasing while for m > 2 it exhibits a 
unique minimum. In between, at m = 2, the function is strictly decreasing with a zero slope at its right 
hand limit. The right panel of the same figure shows how the general shape and the minimal value of 
the function change when m is further increased. In the following paragraphs, we describe the dynamics 
of the system for these three cases of m and illustrate through numerical computations of representative 
examples how this change in the shape oi K translates into the properties of equation (|l.ip . 

5.1. Coherence transition for 1 < m < 2. When 1 < m < 2 the function K = K{c) is decreasing 
on ] — 1, +oo[ (see illustration figure [2] panel (a) for K = 1.5), and the function xi = xi(c) is decreasing 
on ] — 1, +oo[. The bifurcation diagram xi = xi{K) is shown in figure [3] panel (a). 

When K < Kc = 2 the unique equilibrium of (II. ip is the uncoherent steady state (corresponding 
to xi = 0), which is locally and globally stable (label (Pi) figure [3] panel (a)). At K = Kc = 2 a 
supercritical pitchfork bifurcation occurs (see proposition 14. 4p . When K > 2^ the uncoherent steady 
state 2^ is locally unstable (label (P2) figure [3] panel (a)) and equation (jl.ip has a coherent equilibrium 
with xi > 0. When K < K{c = 1) ~ 2.4 this coherent equilibrium is positive on [— vr, tt] (label (Sp) 
figure [3] panel (a)), whereas when K > K(c = 1) this coherent equilibrium is localized: its support is 
strictly included in ] — 7r,7r[ (label (Si) figure [3] panel (a)). 

When 1 < m < 2, solutions u{t,9) of (jl.ip are smooth for all times t > 0. When K < Kc, 
all solutions converge to the uncoherent equilibrium ^ . Figure [5] panels (d) and (e) show one such 
converging solution for K = 1.5. 

When K > K{c = 1) most solutions converge to the coherent equilibrium (Si) of (jl.ip with localized 
support (see illustration in figure [J] panel (f) for K = 6.0). 

5.2. Coherence transition when m > 2. When m > 2 the function xi = xi(c) is decreasing on 

] — 1, +oo[, and the function K = K{c) is decreasing on ] — 1, c* [, it has a unique minimum at c = c* < 1 

and it is increasing on ]c*, +cxd[ (in figure [2] see the lowest curve in panel (a) for m = 3, and panel (b) for 

m = 3, m = 5, m = 8, and m = 15). The value min Km{c), that depends on m, defines the border of 

ce]-i,i[ 

regions (U) and (B) in figure [T) at this value a fold bifurcation gives birth to two (localized) coherent 
equilibria. The equilibria given by (j4.6p with c < c* are coherent, localized and (locally) stable. The 
equilibria given (j4.6p with c < c* are coherent and unstable, for c < 1 they are localized and for c > 1 
they are positive on S. 

When K < min Km{c) the unique equilibrium of (II. ip is the uncoherent steady state which 
ce]-i,i[ 

is locally and globally stable (label (Pi) figure [5] panel (a)). When min Km{c) < K < 2, there is 

C6]-1,1[ 

one coherent localized stable equilibrium denoted by label (Si) in figure [5] panel (a), and illustrated for 
K = 1.84 by the lowest curve in figure [5] panel (d). For the same values min Km{c) < K < 2 there is 

also and one coherent unstable equilibrium that is localized for min Kmic) < K < K{c = 1) (label 

CG]-1,1[ 

(Ui) in panel figure [5] (a)), and positive for K{c = 1) < K < 2. (label (Up) in figure [5] panel (a)). Panel 



11 



(a) Bifurcation diagram for m=1 .5 




c) stable and localized equilibria (S,) when m=1.5 




^ 0.6 
0.4 
0.2 




K = 2.5- 
K = 6.0- 
K = 13.5 




Figure 3. Equilibria bifurcation diagram and equilibria of (II. ip in case m = 1.5. Panel 
(a) show the number of equilibria of p.ip for < < +oo. We see in particular the supercritical 
pitchfork bifurcation occurring at when K — 2. The branch of stable equilibria that appear 
satisfies K — Kc ^ —{m — 2)x\ > in a neighborhood of the bifurcation (see proposition 14. 4|) . 
In panels (b) and (c) the horizontal line represents ^ . Panel (b) shows the shape of coherent 
(stable) equilibria (S'p), for various values of K in ]2,K{c — These equilibria emerge from 
^ when K = 2 and they remains positive until K = K{c — 1). Panel (c) shows the shape of 
localized coherent equilibria {Si), for various values oi K > K{c = 1). When K = K{c — 1) 
the equilibrium is non negative and zero exactly at — n, and the equilibria are localized when 
K > K{c = 1). They converge to a Dirac mass centered at 6* = when K +oo. (see 
proposition (|4.2I ) for analytical formulas for these equilibria). (Details on numerical methods 
used here can be found in the appendix.) 



(c) in figure [5] shows the shape of such localized equilibrium when K = 1.84, and positive equilibria 
are shown for various K{c = 1)K < 2 in panel (b). At K = = 2 a subcritical pitchfork bifurcation 
occurs (see proposition 14. 4p . When K > 2, the uncoherent steady state ^ is locally unstable (label 
(P2) figure [5] panel (a)) and equation (jl.ip has a unique localized coherent equilibrium with xi > 
(label (Si) in figure [5] panel (a)). Two such equilibria are shown in figure [5] panel (d) for K = 2.4 and 
K = 4.7, and in panel (e) for K = 7, K = 17, K = 100 and K = 400. 

When K < K{c*), all solutions converge to the uncoherent equilibrium ^ (region (U) in figure 
[U label (Pi) in figure [5] panel (a)). In the multistability region (B) (see figure [1]), corresponding to 
K{c*) < K < 2, solutions typically converge either to the uncoherent equilibrium (Pi) or to a coherent 
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(d) Convergence to the uncoherent equilibrium 



(e) Convergence to the uncoherent equilibrium 
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(f) Convergence to a localized coherent equilibrium 
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0-6 V ™ = 1.5 
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Figure 4. Convergence of solutions in case m = 1.5. When K < 2, solutions converge to 
the uncoherent equilibrium (Pi). This is illustrated in panels (d) and (e) where we used for initial 
condition "o = ^ + ^ir^ cos(-) with a;i(0) = 0.2, and we have < lO'^ for all t > 10.0. 

Panel (f) illustrates convergence to the localized coherent equilibria (Si) for K > K{c = 1). We 
used K = 6.0, the solution remains smooth and positive before converging, and using the same 
initial condition that in panel (d) it has numerically converged for t > 5. (Details on numerical 
methods used here can be found in the appendix.) 



localized equilibrium {Si). This is illustrated in figure [S] panels (f) and (g). For the same value of 
K = 1.86 and two different initial conditions, we observe convergence to a coherent equilibrium in panel 
(f) and convergence to the uncoherent equilibrium in panel (g). When K > 2 (region (C) in figure [1]), 
most solutions converge to the localized and coherent equilibrium (label (Si) in panel (a) of figure [5]). 
Panel (h) of figure [6] shows such a solution for K = 7. 

5.3. Coherence transition when m = 2. When m = 2 the function K = K{c) is decreasing on 
] — 1, 1[ and constant on [1, +oo[ (proposition 14. 3| see illustations in figure [2] panel (a)), and the function 
xi = xi{c) is decreasing on ] — 1, +oo[. The bifurcation diagram xi = xi{K) is shown in figure [7] panel 

When K < Kc = 2 the unique equilibrium of (jl.ip is the uncoherent steady state ^ (corresponding 
to xi = 0), which is locally and globally stable (label (Pi) figure [7] panel (a)). At K = Kc = 2 a 
degenerate pitchfork bifurcation occurs (proposition 14. 3p . Any u given by n(-) = ^ + cos(-) with 
< xi < ^ is an equilibrium of (jl.ip (label (Np) figure [3] panel (a)). The uncoherent equilibrium is 
recovered when xi = 0, for any < xi < 0.5 the equilibrium is positive on [— 7r,7r], and for xi = 0.5 
the equilibrium is non negative on [— 7r,7r]. When m = 2 the transition from uncoherence to coherence 
is sharp. The family of equilibria connecting to the localized coherent equilibria appear all-at-once 
when K = 2. In particular we have K — Kc = o(x") for any n > 2 in some neighborhood of the 



13 



(a) Bifurcation diagram for m=3 
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(c) Unstable and localized equilibria (U,) when m=3.0 
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(d) Stable and localized equilibria (S|) when m=3.0 
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(e) Stable and localized equilibria (S,) when m=3.0 
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Figure 5. Equilibria bifurcation diagram, case m = 3. Panel (a) shows the number of 
equilibria of (jl.ip for Q < K < +co. We see in particular the subcritical pitchfork bifurcation 
occurring at ^ when K ~ 2. The branch of unstable equilibria that disappear satisfies K—Kc ^ 
— (to — 2)x^ > in a neighborhood of the bifurcation (see proposition 14. 4p . In panels (b) to (e), 
the horizontal line represents Panel (b) shows the shape of unstable coherent equilibria {Up), 
for various values of K in ]K{c = 1), 2[. These equilibria are non negative when K = K{c = 1), 
and positive when ]K{c = 1),2[. They converge to ^ when K ^ 2. The shape of localized 
and unstable coherent equilibria (t/;), that exist for K in ] mini4r,„(c), iir(c = 1)[ is illustrated 
in panel (c) for K = 1.84. These equilibria appear by fold bifurcation at if = m.va.Km{c), and 
they become positive when K > K{c = 1). Panels (d) and (e) show localized coherent (locally) 
stable equilibria (5*;), for various values oi K > min ii'm(c) (see proposition (|4.2I ) for analytical 
formulas for coherent equilibria). These equilibria appear at K = min Km{c) by fold bifurcation, 
they are localized and stable for all K > min Km{c) , and they converge to a Dirac mass a.t 9 ~ 
when K — )• +oo. (Details on numerical methods used here can be found in the appendix.) 
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(f) Bistability : Convergence to a localized coherent equilibrium 



(g) Bistability : Convergence to the uncoherent equilibrium 
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(h) Convergence to a localized coherent equilibrium 
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Figure 6. Convergence of solutions to equilibria, case m = 3. When K < min^„j(c), 
solutions converge to the uncoherent equilibrium (region (U) in figure [Ij. Panels (f) and (g) 
illustrate the bistability phenomena in region (B) of figure [T] (for K in ] min K,yi{c), 2[). Using 
uq = + ^ cos(-) with xi = 0.5, panel (f) shows convergence of one solution to the localized 
coherent equilibria (Si) when K = f .86, and for the same value of K panel (g) shows that with 
different initial data (xi — 0.01) the solution converges to the uncoherent equilibria. When 
K > 2 solutions remain positive and converge to the localized coherent equilibria (Si) (region 
(Ci) in figure[T]), this is illustrated in panel (h) for K = 7. In panels (f) and (h) the solutions have 
numerically converged at t = 15.0 and remain the same for all t > 15.0. In panel (g) convergence 
is much slower, using xi = 0.01 for the initial condition, the solution is still approaching ^ when 
t — 300.0. (Details on numerical methods used here can be found in the appendix.) 



bifurcation point. As soon as K > 2, the uncoherent steady state ^ is locally unstable (label (P2) 
figure [3] panel (a)) and equation (II. ip has a unique localized coherent equilibrium with xi > ^ (label 
{Si) figure [7] panel (a)). 

When K < Kc, all solutions converge to the uncoherent equilibrium When K = 2 the phase 
space is foliated by the stable manifolds of the equilibria «(•) = + cos(-) where < xi < ^, that 
is to say: any solution converge to one of those equilibria, which depends on the initial data uq. Panels 
(d) and (e) in figure [8] show two solutions of (jl.ip for the same value of = 2 and two different initial 
conditions that converge to two different equilibria of the family (Np). When K > 2, solutions typically 
converge to the unique (localized) coherent equilibrium (Si) of (ll.ip . In figure [HI panel (f) show this 
convergence phenomena for one solution when K = 6.0. 

Bifurcation diagrams for m < 2, m = 2 and m > 2 are compared in figure [9l The value of K at 
which stable and localized equilibria (Si) appear depends smoothly on m, it equals m.mKm{c) when 

c 

m < 2, 2 when m = 2 and Km{c = 1) when m > 2. The unstable equilibria (Up) and (Ui) that exist for 
m < 2 converge to the family of equilibria (Np) at m = 2 and K = 2. Similarly the stable and positive 
equilibria (Sp) that exist when m > 2 converge to the family of equilibria {Np) at m = 2 and K = 2. 
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(a) Bifurcation diagram for m=2 
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(c) stable and localized equilibria (S,) when m=2 
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Figure 7. Equilibria bifurcation diagram, equilibria of (II. ip and convergence of solu- 
tions in case m = 2. Panel (a) show the number of equilibria of for < -ftT < +oo. 
We see in particular the degenerate pitchfork bifurcation occurring at ^ when K ^ 2. The 
branch of equilibria that appear at Kc is flat: there is a 1-parameter family of equilibria 
u(-) — ^ + cos(-) where < xi < ^, connecting ^ (for xi = 0) to the onset of local- 
ized equilibria (for xi = ^) (see proposition I4.3p . Positive coherent equilibria of his family (Np) 
are shown in panel (b), for K = 2 and various values of xi = Xi{c), and confirm that the 
equilibrium u converge to ^ has c — )■ +oo and xi — ?> 0. Panel (c) shows the shape of coherent 
equilibria {Si), for various different values of K > K{c — 1) (see proposition (|4.2I ) for analytical 
formulas for these equilibria). These equilibria are smooth and localized for &\\ K > 2, and they 
converge to a Dirac mass aX — Q when K — > -|-oo. (In panels (b) and (c), the horizontal line 
represents details on numerical methods used here can be found in the appendix.) 



5.4. Coherence transition when m — )• -|-oo. The limit m — )• -|-oo is a limit of slow diffusion in 
(jl.ip . The bifurcation scenari in the limit m large is illustrated in figure 15.41 The fold bifurcation point 
K = mine -^(c) coonverge to the pitchfork bifurcation point K = 2, and the interval K G [mine K[c), 2[ 
in which equation (ll.ip is bistable vanishes in the large m limit, while the coherent equilibria branch 
{Si) converges pointwisely to the unstable uncoherent equilibrium branch {P2) (see illustrations of these 
convergence phenomena in figure . 
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(d) Convergence to a positive colierent equilibrium 
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(e) Convergence to a positive colierent equilibrium 
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(f) Convergence to a localized coherent equilibrium 
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Figure 8. Convergence of solutions in case m = 2. Panels (d) and (e) show convergence of 
solutions to different coherent equilibria when K — 2. We used uo{9) = 57 + 7 cos(6') + ^ cos{9) 
with a = 0.3, b = 0.1 in (d), and a = 0.4, b = 0.2 in (e). The solutions converge to m(-) = 
^ + ia;iC0s(-) with xi ~ 0.3 in (d) and xi ~ 0.36 in (e). When K > 2 solutions remains 
positive and converge to the localized coherent equilibria (Si) of (|l.f p . This is illustrated in 
panel (f) where the solution has numerically converged at t = 5.0 and remains the same for all 
t > 5.0. (In panels (d) to (f), the horizontal line represents 
used here can be found in the appendix.) 



details on numerical methods 



6. Discussion 

In this work, we have shown that nonlinear diffusion alters the dynamics of eq (jl.ip in two key aspects. 
One relates to the very shape of coherent equilibria that become localized when advection is strong. The 
other modification is that for m > 2, the transition to coherence scenario is through a subcritical rather 
than a supercritical pitchfork bifurcation, thereby leading to existence of a multistable regime that has 
no counter part in the case of linear diffusion. Central to this impact of nonlinear diffusion is the fact 
that the pitchfork bifurcation for m = 2 is highly degenerate and acts as an organizing center. In the 
following paragraphs we first discuss the implications of this organizing center and then the generality 
of our results. 

Schematically, organizing centers are degenerate bifurcation points where a number of qualitatively 
distinct regimes come to meet. For eq (|l.ip the highly degenerate pitchfork bifurcation taking place 
at {K, m) = (2, 2) plays this pivotal role. We claim that small perturbations of eq (11. ip at this point 
can produce arbitrary one-dimensional dynamics of the interval. Indeed, a full continuous interval of 
equilibria exists at this point. This fact combined with the claim (based upon numerical investigations) 
that the linearized operator at these equilibria possesses a spectral gap suggest that this equilibria 
interval is a normally hyperbolic invariant manifold. This property in turn implies that this manifold 
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Bifurcations diagrams for m<2, m=2, m>2 
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Figure 9. Transition in bifurcations diagrams at m = 2. Tlie value of K at which stable 
and localized equilibria {Si) appear depends smoothly on m, it equals miTiKm{c) when m < 2, 

c 

it equals 2 when m = 2 and i^m(c = 1) when m > 2. The unstable equilibria {Up) and ([/;) that 
exist for m < 2 converge to the family of equilibria {Np) at m = 2 and K = 2. Similarly the 
stable and positive equilibria (5*^) that exist when m > 2 converge to the family of equilibria 
{Np) at m = 2 and ^ = 2. 

Bifurcations diagrams for m large 
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Figure 10. Bifurcations diagrams for large m. In the large m limit, the fold bifurcation 
value K = minc-fi'(c) converges to K = 2, and the uncoherent equilibria {Ui), {Up) exist in a 
vanishing intervall of K (corresponding to the bistability region {B) is figure [T]). The coherent 
localized stable equilibria branch {Si) converge pointwisely to the uncoherent equilibria branch 
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would persist under small perturbations of eq (jl.ip . Finally, the fact that dynamics on the equilibria 
set is trivial allows one to construct ad- hoc small perturbations to reproduce arbitrary dynamics. 

A related construction exists for linear diffusion where it has been proved that small perturbations 
can produce arbitrary dynamics of the circle |21) . While the main arguments, i.e. the existence of a 
normally hyperbolic contiuum of equilibria are similar in both cases, there exists a main difference. For 
linear diffusion, the continuum of equilibria does not exist if one considers even solutions, whereas for 
m = 2, it does so. In other words, if we relax the constraint on solutions beeing even, we can expect 
small perturbations of eq (jl.ip to produce arbitrary dynamis on the disk (and not only the interval or 
the circle). In this sense, our results suggest that the dynamics with nonlinear diffusion are richer than 
those with linear diffusion and give a precise meaning to this property. 

In our work, results were presented for the case where the nonlocal advection term takes on the 
form of a convolution with a sine function. Here, we discuss extentions of our results to other forms of 
advection. 

For the sake of simplicity, we have restricted our analysis to the coupling term J * u = /g sin(- — 
ip)u{(p)df. Results about convergence of equilibria can be extended to more general coupling terms, 
including the Maier-Saupe potential J *u = jgSin2(- — (p)u{(p)dLp for example and more general non 
local convolution with an odd function. In the case of Maier Saupe potential the bifurcation scenario 
remains very similar to figure [TJ The main difference is that coherent equilibria have two maximum on 
[— TT, vr], and localized coherent equilibria have non-connected support of the form Jvr — e, tt + e[U] — e, e[. 

Polynomial advecting terms have been considered in the context of Keller-Segel model with degenerate 
diffusion for example see |39l l38l HQ] . The corresponding equations 

dtu = a|(n™) + de {n't J *u), t > 0, G [0, 27r], 

u{Q,e) =uoie) e€[0,2Tr], 

u(t,0) =n(t,27r) t>0, 

d0u{t,O) =deu{t,2TT) t>0, 

with 1 < q < m can be treated as the case q = 1 and the bifurcation scenari are the same, with value 
m changed to m — {q — 1). Here we simply show that proposition 14.21 on the equilibria set of (jl.ip can 
be directly extended to the following corollary 

Corollary. Let u be a stationary solution of ( fg. 1\) with 1 < q < m. Then we have 

m 



u^-i{e) = ——LKxi [cos(0) + c]. 



where xi and K are given by 



Jm-(q-l) m 1 1 m 

xi = - — ^ K - - 



' rn-qx^ C-l-l) " ^ (<?-!) C-Vd 

for c g] — 1, +oo[, where Im and Jm are defined as in proposition^^ 
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Appendix 

Figure [3] panels (a) to (c), figure [7] panels (a) to (c), and figure [5] panels (a) to (e), were made using 
formulas (|4.4p and (|4.6p that give xi, K and the coherent equilibria u as functions of c. The functions 
Im{c) and Jmic) were approximated by the classical trapezoidal rule with space discretization Ag = 1Q~^ 
for c g] — 1, 2[ with discretization Ac = 10~^. The coherent equilibria are stable exactly when c < 1 and 
K > min i^(c). The coherent equilibria are localized exactly when c < 1. The bifurcation diagram in a 
neighborhood of K = 2 and xi = is obtained in the limit c — )• -|-oo. Since we have proven analytically 
that K — 2 = 0{x1) in that limit, the have extended the curves obtained with c g] — 1, 2[ by straight 
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lines for c > 2. Meanwhile the values min K(c) were recorded for each m, and then used to draw the 

CG]-1,1[ 

border between regions (U) and (B) in figure [TJ In the same figure the borders between regions (Bi) 
and (Bp) and regions (C;) and (Cp) were obtained for K{c = 1). 

Simulations of solutions of (|l.ip shown in figures in the main text have been done using a finite 
difference scheme to discretize the diffusion operator 



dlu^iO) 
e advectic 
d${uJ * u) ^ Kxi- 



a symmetric scheme for the advection term 

u{9 + Ag) sm{9 + Ae) - u{e - Ag) sm{9 - Ag) 



2Ae 

and a explicit Euler scheme for the time derivative. Unless something else is explicitly mentioned, 
we have used the initial condition uq{0) = ^ + ^cos(0). We mention that the scheme slightly 
unstabilize the uncoherent equilibria even when K < Kc- It also increases diffusion phenomena and 
induces some numerical instabilities at the border of the support of solutions when this support is 
strictly included in ]0,27r[. Numerical results have been checked using (Ag, At) = (2.0 10~^, 10~^) and 
{Ag, At) = (2.0 10~^, 10~^), convergence has been checked when Ag, At — )• 0, and the convergence rate 

is 0(A,) + 0(At) + 0((^). 
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